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Abstract 

Perturbative QCD predicts that at sufficiently large Q 2 and small x nucleon structure functions 
should exhibit scaling in the two variables \ lnlnQ 2 and y^ln ^/ln InQ 2 , provided only that 
the small- a; behaviour of the input to the perturbative QCD evolution is sufficiently soft. We 
derive these asymptotic results by writing the gluonic Altarelli-Parisi equation at small x as 
a two-dimensional wave equation, which propagates the gluon distribution from its boundaries 
into the asymptotic region. We then show that the existing experimental data on Ff^x, Q 2 ) 
from HERA provide a remarkable confirmation of both of these scaling predictions. The 'hard' 
pomeron, which does not scale, is thereby excluded by more than three standard deviations; only 
a very small admixture of it is permitted by the data. We propose that existing and future data 
from HERA should be binned in the two scaling variables, in order to facilitate the search for 
small scaling violations. 
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The main evidence for perturbative QCD in deep inelastic scattering at large x comes 
from the observation of small violations of Bjorken scaling [|I] in the structure function 
F2(x;t) in the Bjorken limit of large t = lnQ 2 /A 2 at fixed x. These scaling violations 
may be accurately predicted using the operator product expansion or equivalently the 
Altarelli-Parisi equations . In this letter we will show that in the double limit of large 
t and small x perturbative QCD also makes definite scaling predictions, provided only 
that the gluon distribution at scales of order 1 GeV is reasonably soft at small x (i.e. 
that the pomeron intercept is close to one). The appropriate scaling variables are the 
geometric mean a ~ ^/ln ^ In t and the ratio p ~ \J In ^ / In t ; the gluon distribution then 
scales asymptotically in both a and p in the double limit of large a at fixed p and large 
p at fixed a. These simple predictions, which are asymptotically parameter free up to 
an overall normalization factor, may be directly tested by comparison with experiment. 
Indeed, we will show that recent data from HERA HJ| are in remarkable agreement 
with this asymptotic behaviour. Furthermore, since the double asymptotic scaling would 
necessarily be spoiled if the pomeron were hard, we will be able to use it as a means to 
discriminate experimentally between soft and hard pomerons. 

Before actually proceeding to the comparison of double asymptotic scaling behaviour 
with the data, we first derive it in a simple but, we believe, instructive way by reducing 
the singlet Altarelli-Parisi equations at small- x and large t to a wave-like equation for 
the gluon distribution, which may be solved exactly. We then derive the asymptotic form 
of this solution for large values of the scaling variables a and p, and different classes of 
boundary conditions, namely 'soft' and 'hard'. We explain in what sense solutions with 
soft boundary conditions exhibit double asymptotic scaling, establish the kinematic region 
in which this scaling behaviour sets in, and show the way in which it is corrupted if the 
boundary conditions are 'hard'. We derive a simple relation between F2 and the gluon 
distribution in the asymptotic region, and finally proceed to a direct comparison with the 
experimental data. 

We will be concerned specifically with the singlet contribution to the proton structure 
function 

F*(x;t) = &q(x;t) + F^ s (x;t) 

Uf (!) 
q(x;t) = (ftfo^ + ftte*)), 

i=l 

where rif is the number of active flavours. At small enough x the nonsinglet contribution 
F^ s (x;t) is negligible and can be ignored. At small x and large t the singlet quark 
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distribution q(x; t) is essentially driven by the generic instability of the gluon distribution 
g(x; t); determining the latter thus allows us to compute F% in this kinematic region. 

To see how this works, consider the singlet Altarelli-Parisi equations ||, which de- 
scribe perturbative evolution of g(x; t) and q(x; t) in the kinematical region displayed in 
fig- B 

3 g(x; t) = ^ t ^[P g9 (^g(y; t) + P gq {fMy; *)], 



dt 2tt , , ., 

" x u (2) 

!L~f- as ^ f dy > 

dt 



q(x;t) = ^- [ d ^[2n f P qg ^)g{y-t) + P qq ^) q {y-t)l 

z7r Jx y 



where a s (t) = J^- + O(-p), (3q = 11 — |n/. Solving these equations sums up all leading 
logarithms of the form a™ ln m Q 2 ln n - with n < m; terms with fewer powers of In Q 2 should 
be negligible if t is large enough for perturbation theory to be trusted (t > to say, where 
a s (t ) < 1). 

Now evolution at small x (x < xq , say) will be dominated by the lowest moments of the 
splitting functions; this may be demonstrated explicitly by taking the Mellin transform of 
the evolution equations (||) , solving them, and then using a saddle-point evaluation of the 
inverse Mellin transform of the solution |||7[]. In practice, this is equivalent to expanding 
the Mellin transform of the one loop splitting functions about their largest singularity in 
moment space (which is when the moment variable n — > 0) and discarding terms which 
vanish as n — > 0[j^,^l. Inverting the transform again, this leads to the approximate form of 
the splitting functions 

p 9 M + ^W 1 - *), PM « h - 25(1 - z), 

P qg (z) « |5(1 - z), P qq (z)^0. 

When substituted into @, the terms in (H) which go as z~ x will sum all the double 
logarithms, while the terms proportional to 5(1 — z) will approximately sum the remaining 
single logarithms^; in terms of ladder diagrams, the former terms account for rungs which 
are strongly ordered in both x and t, the latter for rungs which are only strongly ordered 
in t. If only the former terms were retained one would get the so-called double-leading 
logarithm approximation. 



1 This approximation could be refined by keeping higher powers of n in the Mellin transform 
of the splitting function. This leads to further terms in @ involving derivatives of the delta 
function. 
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Defining for convenience two new variables 



(so that at a; = xo ^ = 0, and at i = to C = 0), and writing 

G(^C) = xg(x]t), Q(£,C) = xq(x;t) = ^ x (q^x; t) +qi(x;t)), (5) 

i 

then substituting (||) into (||), and differentiating the first equation with respect to £, the 
singlet evolution equations assume the simple form 

() " o + o - 7 2 g(^, o = o - o, (6) 



^Q(e,o = ^G(e,c), (7) 

where 7 = ^/12//3q and <5' = (11 + ^-)/(3q. These equations should be valid throughout 
the first quadrant of the (£, C) plane ( fig. ||), outside the multiple rescattering region 
£ 3> e 2 ^ [Q, up to corrections of order xoe - ^ and a s (t )e~^ .i 

The inhomogeneous term on the right hand side of (^|) is negligible due to its small 
coefficient and may be dropped. However, this is not necessary: if the matrix of splitting 
functions is diagonalized before the Mellin inversion,! this term is in effect absorbed into 
a small shift in 8', which becomes 5 = (11 + -^f) / Po- The gluon evolution equation (|]) 
then becomes homogeneous, 

°~ +5|:-7 2 teO = 0, (8) 



while the quark evolution equation (|7|) remains unchanged (up to subleading corrections) . 



2 Corrections to the expansions (||) of the one loop splitting functions will lead to contributions 
to the parton distributions which are suppressed by powers of p~ x (where p is as defined in (11) 
below). Furthermore it may be shown by explicit computation that the two loop corrections to 
(|3|) have little effect on the final form of F2 provided a <^ 1 and 1 & p & 3. However for n & 0.35, 
which in practice 3, there are indications that perturbation theory may begin to break 
down |0. 

3 Although in general the diagonalization introduces an extra singularity (a branch cut) this 
is not a problem here since we expand about n = 0. 
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We next make the essentially trivial observation that (^|) is in fact a two dimensional 
wave equation, with 'time' i(£ + £) and 'space' i(£ — £)). This allows us to immediately 
infer many useful properties of its solutions, and thus deduce directly many well-known 
qualitative features of parton evolution. In particular: 

(i) The equation is essentially symmetrical in £ and £, so the gluon distribution xg(x; t) 
evolves ('propagates') equally in both x and t (up to the (small) asymmetry induced 
by the 'damping' term proportional to 5). Any further asymmetry in £ and £ must 
thus come from the boundary conditions. 

(ii) The propagation is 'timelike', into the forward 'light-cone' at the origin (£, Q = (0, 0), 
along the 'characteristics' £ = constant and ( = constant (see fig. |J). 

(iii) At a given point (£, (), G(£X) depends only on the boundary conditions contained 
within the backward light cone formed by the two characteristics through (£, Q fig. 0. 
This means that to calculate xg(x;t) it is unnecessary to know what happens at 
smaller values of x or larger values of t. 

(iv) Because the equation is linear, contributions to G(£, Q from different parts of the 
boundaries are simply added together. 

(v) Since the 'mass' term is negative (7 2 > 0), the propagation is 'tachyonic'; this means 
that G(£, is unstable, growing exponentially rather than oscillating. 

(vi) Since 5 > the damping term ensures that, at fixed £ G(£, () eventually falls with 
increasing (. 

Setting G(£, () along two characteristics (for example £ = and ( — 0) gives a well- 
posed problem (the 'characteristic Goursat problem' |TJ]) provided only that the boundary 
conditions are compatible at the point where the boundaries meet. Since the Green's 
function for the wave operator (|8p may be expressed in terms of the Bessel function 

00 (±z 2 ) n 1 
\n\) \/2nz 

the general solution to this problem is simply 

< , a 



G(e,C) = Io(27VeC)e- 5C G(0,0) + J d$%(?yy/(Z - Z>)C) e - s <-^G{?, 0) 

+ J* dC% (2 7V / e(C-C / ))e 5(C '- C) (^(0, C') + SG(0, (')) . 



(10) 



The asymptotic behaviour of G away from the boundaries will depend in general on 
whether the boundary conditions are compact or noncompact, i.e., on whether the starting 
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distributions are soft or hard at large £ on the £ = boundary and large £ on the £ = 
boundary. Due to the linearity property, we can consider each boundary separately, and 
then add the various contributions. We consider first the simplest case in which both 
boundary conditions are soft, so g(x;t Q ) ~ x" 1 as x — > (i.e. G(£,0) ~ const, as £ — > oo), 



the traditional behaviour inferred from the intercept of the 'soft' pomeron ||12|| , while 
g{xQ\t) ~ t~ s as t — > oo (i.e. G(0, C) ~ e _<5 ^ as £ - ► oo). Since in this case 'waves' are 
only generated close to the origin, far from the origin we may determine the behaviour of 
the solution fllCf ) by means of a multipole expansion, i.e., expanding the Bessel functions 
under the integrals in eq. ( p!Q[ ) in powers of and ^- respectively. Defining the new scaling 
variables 




and considering the limit o — > oo at fixed (but reasonably large) p (see fig. |]),@ we then 
find asymptotically the generic behaviour 

G(a,p) ~ / (27a)e- 5(fj/p) (G'(0,0) 

(T — + OO I 

"fM^s «p {^-*(5)}(i + o(i)). (12) 

where y = 7/p, and in the second line we used the asymptotic behaviour (|^) of the 
Bessel function. The normalization factor Af is just that of the boundary condition!, while 
f g (0) = 1, but f g (y) will depend on the particular details of the boundary condition, 
while being generally rather smoothJl The asymptotic behaviour (|l~2| ) is thus, as expected, 
independent of the detailed form of boundary condition, up to corrections of order 1/p, 



4 More precisely, we consider the limit a — > 00 along any curve such that also p — > 00, such as 
for example the curve £ oc C 1+e with e > 0. 

5 The boundary condition g(x;to) will be normalized so that the second moment (the gluon 
momentum fraction) is unity; the boundary condition g(xo,t) is then normalized by continuity 
at (xo,to). When comparing to the data, the overall normalization will thus give the gluon 
momentum fraction. 

6 For example if xg(x;t ) = Af(l - x) , x g(xo;t) = Af(l - x ) /3 (^) 5 , then f g {y) = (1 - 
*of + Pxo ELI J ( S:l- n y. \7:Tv - l -^TT- y + °(*°) as x - 0. 
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which are anyway of the same order as corrections due to the terms lost by the truncation 
([D of the expansion of the splitting functions. 

If instead the boundary conditions are noncompact, it is not difficult to see from (|T0| ) 
that the form of G(£, on the boundary will tend to persist into the interior of the (£, 
plane; the asymptotic behaviour close to the boundary will then not be universal, but 
depend on the particular boundary conditions adopted. Consider for example xg(x; t ) ~ 
J\f'x~ x , (phenomenologically A ~ 0.08 [0). Then if p > ^, the first boundary integral in 
eq.(|lOD is dominated by a nontrivial saddle-point, and gives the asymptotic behaviour 

G(a,p) ~ A/-'exp{Aap+(^-5)(^)}(l + 0(i)), (13) 

a — * oo " 

to be added to (0). Very close to the boundary (p <i 15) the power-like growth at small x 
is thus preserved by the evolution to t > to- However further away from the boundary, in 
the region p ? the more gentle growth ([TJD is unchanged, modulo a slight increase in J\f. 
Similarly, if the boundary condition on the other boundary is hard (say G(0, Q ~ e~ s ** 
as C — > oo, with 5' < 5), this too will propagate into the region p ^ ^ , in effect 
reducing the value of 5 close to the boundary. Again away from the boundary the universal 
behaviour ([12]) is unaffected. 

Thus with only the rather conservative assumption of reasonably soft boundary con- 
ditions (i.e. for sufficiently small A and 5') perturbative QCD predicts a universal growth 
( |12|) in the gluon structure function at large t and small x, faster than any power of In 
1 nit slower than any inverse power of x. This means that away from the boundaries pertur- 
bative QCD allows us to predict in a parameter-free way (i.e. independently of the precise 
form of the boundary conditions) the asymptotic behaviour of the gluon distribution, and 
hence (as we will see below) of the F2 structure function, throughout a very large region of 
the (x,t) plane fig. |l|. Although it has been known now for twenty years 0], hitherto this 
remarkable consequence of the perturbative QCD evolution equations seems not to have 
been taken very seriously. Here we hope to remedy this. 

We begin by noting that the universal asymptotic behaviour (|12D actually implies scal- 
ing in both a and p: lnG(a, p) at fixed p is asymptotically a linear function of a, with slope 
independent of p (up to terms which vanish as 1/p), while at fixed a it is asymptotically 
a flat function of p. In other words a" 1 In G(cr, p) is asymptotically independent of both a 
and p, up to corrections which vanish as 1/cr and 1/p. This 'double asymptotic scaling' 
is a direct result of two generic properties of wave propagation from compact boundaries: 
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at large distances from the boundaries, the details of the boundary conditions are washed 
out (a scaling)@, and, apart from the damping term (which is asymptotically subleading) 
the propagation is isotropic (p scaling). 

This double scaling property is peculiar to the universal behaviour (fT^ ) derived from 
reasonably soft boundary conditions; instead, if one of the boundary conditions is very 
hard, its influence will spread into a large part of the (£, Q plane, spoiling both p and a 
scaling in this region in a way that will depend explicitly on the particular form of the 
hard boundary condition. For example consider the Lipatov 'hard' pomeron ||15|| , which 
grows as a power of - as x — > at fixed t: xg(x;to) ~ J\f'x~ x , with A 12 ^ n2 a s (t ) « ~. 
This behaviour is supposed to be valid at some intermediate t: it is obtained by summing 
leading logarithms of the form ck™ ln n ~, terms containing powers of InQ 2 being ignored; 
t must be small, but not so small that a s (t) is too large. It will thus serve as a boundary 
condition for Altarelli-Parisi evolution into the region t > to. It will then lead to power- 
like growth of the form (|13|) throughout the region p ~ Q~ , which must be added to the 
weaker growth (|T2p over which it will dominate at large p. Although the structure function 
still grows at small x, it now increases so strongly that double scaling is spoiled. 

In what remains of this letter we will compare the double asymptotic scaling pre- 
diction (|12D with all the available small-x data from the two HERA experiments @||, to 
see whether it is supported empirically, or whether there is instead evidence for other less 
conventional physics (and in particular the 'hard' pomeron) at small-x. Firstly we must es- 
tablish that there are data in the asymptotic region where double scaling might apply. We 
thus begin by examining the approach to asymptopia as the general solution (|10|) evolves 
away from the boundaries. We choose xo = 0.1, which should be sufficiently small for (|^) 
to be a good approximation. The starting scale is chosen as Qo = lGeV: not so small that 
perturbative evolution has broken down, but sufficiently small that the behaviour at small 
x is given by a nonperturbative boundary condition (or, in the case of the hard pomeron, 
sufficiently small that logarithms of Qq may be reasonably ignored). 

In fig. |3] we display a contour plot of G(£, () computed from the solution (|T0| ) with, 
for definitness, the soft boundary condition G(£,0) = jV(1 — xoe~^) 9 compared with the 
solution computed from the hard boundary condition G(£, 0) = Af'x 1 ' 2 e^^ 2 (l — xoe - '') 9 . 



7 The position of the boundaries is also unimportant asymptotically, so even if to is taken to 
be very small the same generic behaviour will be found at small x and large t; as we will see, this 



is sufficient to explain the success of the predictions of ref. [14|. 
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On the £ = boundary we take in both cases G(0, C) = G(0, 0) as (approximately) required 
empirically when xq = 0.1. The plots clearly show that the asymptotic behaviour is 
attained very rapidly in the soft case, around a ~ 1; the damping term and the propagation 
from the lower boundary distort the curves somewhat at small values of p, as expected. 
In the hard case the asymptotic behaviour ([R|) clearly holds only for large 

3, while in 

most of the rest of the region the full solution is a compromise between the two sorts of 
behaviour. In practice, this means that in much of the HERA kinematic range the evolved 
hard solution is actually rather similar to the unevolved boundary condition imposed at 

C = o. 

It is now clear that, because almost all of the HERA data (and even some of the older 
NMC data) lie in the region a G [1,2], p e [1,3], well inside the asymptotic kinematic 
regime of the soft solution, the universal behaviour ([12]) should be easily tested, while the 
hard behaviour should be easily distinguished from it. This is demonstrated explicitly in 
fig. |] where we compare the scaling expected from the soft pomeron boundary condition 
with the hard pomeron prediction (with A = h). Not only does the hard pomeron generally 
lead to more gluons in any given range of x, but the number rises rather more steeply at 
smaller values of x. Thus at fixed p (fig. [|i), even though lnG(a", p) is in both cases an 
(approximately) linear function of a in the hard case its slope gets steeper for larger values 
of p, while in the scaling case it is universal. Correspondingly at fixed a, in the soft case 
G(a, p) becomes flat as p grows, while in the hard case it rises rather rapidly. This is best 
exhibited by first performing a simple rescaling of G(a, p), by a factor 

R G (a, p) = exp(-2 7 (7 + 5{a / p) + \ lna), (14) 

so that the curves with different values of a lie on top of one another (fig. |]ii). Besides 
showing the clear distinction to be made between the predictions derived from soft and 
hard boundary conditions, these plots also confirm that for compact boundary conditions, 
the asymptotic double scaling regime is indeed reached rather rapidly, at moderately small 
values of a and p. 

Before coming finally to a direct comparison of this predicted behaviour to the data 
]4|,[5| we must still relate the calculated gluon distribution to the proton structure function 
F% measured in the experiments. This is actually rather simple at small x and large t, 
since all that is necessary is to integrate (0): 

Fl = C) = 0) + Sf7 2 / dCG(e, O- (15) 

8 



With compact boundary conditions, the asymptotic form of F% may then be deduced 
directly from that of G, by substituting ([12]) into (|T5|), to give (whenever p > 5/7) 



~ ^V/ fl (v)G(<T,p)(l + 0(±)), (16) 

a — ► 00 

where again y = 2, and — * 1 as y — > 0.@ A similar relation holds for the asymptotic 
behaviour (|I3p, but with y = X when p > j. So we can actually use QIfi| ) as a good 
approximation throughout the asymptotic region for any type of boundary condition if we 
set y = ■sj (7/p) 2 + A 2 . Apart from the subasymptotic factor of j/p the asymptotic growth 
of the gluon distribution drives a similar growth of the singlet quark distribution, and thus 
of the structure function F%. 



Using (|T6|) to determine a slightly different rescaling factor 



R F (a,p)= [fej] _1 i2o((T,p). (17) 

we may thus compare the a and p scaling plots fig. ^| to all the available data from ZEUS 
0] and HI ||. The result is shown in fig. |5|. It should be immediately apparent that the 
agreement with the asymptotic behaviour (|T^) derived from the soft pomeron boundary 
condition at to is remarkably good, especially when it is remembered that the asymptotic 
prediction is parameter free, apart from the precise value of the normalization. The hard 
pomeron boundary condition on the other hand gives results which are a long way from the 
data, even when its normalization is reduced. Furthermore, the a plot provides a nontrivial 
test of perturbative QCD, since the fitted slope agrees very well with the computed one. 
Also included in fig. |5|i are the data from the NMC collaboration jn3|; remarkably those 



which fall within the asymptotic region agree very well with the line fitted to the HERA 
data alone, despite their much smaller errors. 

We may assess quantitatively the quality of the fit of the predictions to the HERA 
data by computing the x 2 of the hard and soft structure functions (Ff ) soft and (F$ ) hard 
(obtained using (|16[) from the gluon distributions displayed in fig. ^|a) and b) respectively), 



In fact the function f q may be determined by integrating the second Altarelli-Parisi equation 
dH) with the second term set to zero (since its contribution is subasymptotic) and P qg (z) = 

\ [z 2 + (! - z f] ! this g ives fq(v) = (! - ^zV)" 1 2{i+t)(lVy){z+y) ■ Since we alread y have a similar 
uncertainty from f g in (fi"2|), due both to our small x approximation (||) to the gluon evolution 

equation and ignorance of the boundary conditions, this result will not be used here, though it 

may be useful for extracting the gluon distribution directly from the data. 
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N s 


N h 


x 2 


a) 


0.33 ±0.01 





18 (52) 


b) 





0.14 ±0.01 


114 (282) 


a)+b) 


0.31 ±0.03 


0.01 ±0.01 


17 (50) 


c) 





0.36 ±0.01 


121 (285) 


a)+c) 


0.30 ±0.05 


0.03 ±0.01 


17 (49) 


d) 





0.17 ±0.01 


67 (168) 


a)+d) 


0.30 ±0.05 


0.02 ±0.01 


17 (50) 



Table: The fitted normalizations N s and Nh and the associated x 2 s (those in parenthesis 
using only the statistical errors in the experimental data). The different cases considered are 
a) soft pomeron b) hard pomeron c) hard pomeron (unevolved) d) 'intermediate' pomeron, 
and then the linear combinations a) + b), a) + c) and a) + d). 

with the overall normalization of each left as a free parameter. The results are displayed 
in the table. The \ 2 of the soft pomeron is rather good and fully compatible with the 
experimental errors in the data: there are 38 data points, with presumably some correlation 
in their systematic errors, and our % 2 is indeed rather low if we combine the statistical 
and systematic errors for each point in quadrature, but somewhat high if we use only 
the statistical errors. The fit of the hard pomeron to the data is in contrast very poor, 
excluding it at the level of at least five standard deviations. Notice that the best-fit value of 
the normalization is in fact the fraction of the nucleon's momentum carried by the gluons 
at the starting scale £ = 0; the value found in the case of soft boundary conditions is thus 
in good qualitative agreement with phenomenological expectations, whereas that for the 
hard boundary condition is typically too small by a factor of two. As a cross-check, we 
also fitted the linear combination N^F^ ) soft + N^F^ ) hard , so that both the absolute and 
relative normalizations of the hard and soft components are left as free parameters. As 
may be seen from the table, only a very small admixture (less that about 5%) of the hard 
component is permitted. 

We have also determined the corresponding results for an unevolved hard pomeron 
(again with intercept A = |). This would presumably be appropriate if the Lipatov 
equation were somehow to remain valid even when t ^> to, i.e. in the HERA region.i. It 

9 Amusingly this unevolved hard pomeron may be imitated by adding a negative damping term 
—XdG/dC, to the left hand side of the gluon evolution equation (H), thereby including a™ ln n ^ 
terms in the evolution, and then evolving the soft pomeron from t = to- 



10 



may be seen from the table that in fact this makes very little difference, the improvement 
in fit as compared to the evolved hard pomeron being at best marginal. Since these are 
two extreme scenarios, they bracket the intermediate possibility that the value to at which 
the Lipatov boundary condition should be imposed is somewhere in the middle of the 
HERA region. In the same spirit we also include an 'intermediate' pomeron with intercept 
A = 0.35, since this is supposed to be the minimum value expected from the Lipatov 



equation [IT]. Although this fits a little better than the hard pomeron, it is still clearly 
excluded. 

Finally, we have checked that the x 2 of the soft pomeron solution is not changed 
significantly if the value of to is reduced, or if if the value of A is varied within present 
experimental uncertainty. Of course, if to is raised too much, so that the data are not yet 
in the asymptotic region of the soft solution, then x 2 starts growing and becomes almost 
as bad as that of the unevolved hard solution: while the data clearly exhibit asymptotic 
double scaling, this necessarily implies a growth at small x for fixed t which could easily 
be misinterpreted as evidence for the hard pomeron. 

We conclude that there is simply no room for more than a very small admixture of 
the hard pomeron when x <^ 10 -4 . On the contrary, the data agree surprisingly well with 
the universal double scaling behaviour which only applies in the case of reasonably soft 
boundary conditions. These results seem thus to agree with the recent observation |[L8|| 
that if the integrated kr in the Lipatov equation is cut-off in the ultraviolet, as required 
kinematically, A is very significantly reduced for reasonable values of t , and it becomes 
very difficult to distinguish the hard perturbative pomeron from the soft nonperturbative 
one. 

Since all existing published data is binned in x and Q 2 , it was necessary when drawing 
our scaling plots fig. |5| to use all the data in each of the plots, rescaling with the function 
Rf to ensure that if the data exhibited scaling, they would lie on a single straight line in 
each of the plots. A rather better test of asymptotic double scaling could be made if the 
raw data were binned in p and a: plots of ln-Ff (c, p) at fixed p and a could then be drawn 
to display both p-scaling and a-scaling. Whereas the approach to a flat (p-independent) 
line of the plots at fixed a would indicate that asymptotic scaling has been attained, the 
value of the asymptotic slope of the plots at fixed p would be given by the universal value 



27 = 12/ ^/33 — 6rif/N c , a parameter-free prediction of perturbative QCD which depends 
only on the number of colours and flavours. 
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On the other hand, as much more data becomes available, it should become possible to 
sharpen the tests of double scaling to the point where the subasymptotic scaling violations 
at small a and p can be observed and measured; this would be interesting as some of the 
leading violations (those which are included in (|14|)) are universal, while the subleading 
violations could be calculated quite precisely using the known parton distributions at large 
x and the full evolution equations. It should then be possible to extract (in a rather clean 
way) a s from the HERA data. More excitingly, one could also search systematically for 
scaling violations at large a and p due to new effects such as the hard pomeron or parton 
recombination. The accurate determination of structure functions at small x will then 
not only provide a precise and highly nontrivial test of perturbative QCD, but could also 
be used to search for new physics at small x. 
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10 Indeed the /o-scaling plot fig. |5|b) does seem to indicate a small (but not as yet statistically 
significant) rise for p <^ 3. 
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Figure Captions 

Fig. 1. Region of validity of the Altarelli-Parisi equation in the (x,t) plane. The disal- 
lowed regions are t < to, where perturbation theory breaks down because a s (t) 
becomes too large, and x & x r exp(— a s (to) 2 /a s (t) 2 ), where recombination ef- 
fects may become important @. The small- x evolution equations may be used 
above the dashed line. The dotted region with the ? is that in which the Lipatov 
equation may become applicable. 

Fig. 2. The (£,C) plane, showing the backward light cone at the point (£',C')> curves of 
constant a (the hyperbolae) and lines of constant p. 

Fig. 3. Propagation of G into the £-£ plane; a) soft pomeron, b) hard pomeron. A scatter 
plot of the NMC and HERA data is superimposed for reference: the two types of 
crosses are NMC data [|16| (taken with 90 GeV and 280 GeV beams respectively) , 
the diamonds are ZEUS data f§, and the squares HI data ||. 

Fig. 4. Scaling plots of the function G(a, p) displayed in fig. [|. The solid lines correspond 
to fig. |3|a) (the soft pomeron) and the dashed lines to fig. |3|b) (the hard pomeron) ; 
the dotted lines correspond to an unevolved hard initial condition. 

i) ln(7(cr, p) is plotted versus cr; the lower line of each pair has p = 1.4 and the 
upper one p = 3.2. 

ii) The rescaled function Rg(c, p)G(a, p) is plotted versus p; the lower line of each 
pair has a = 1.2 and the upper one a = 2.1. 

Fig. 5. Comparison between the experimental data and the predictions for ^(a, p) com- 
puted from G(a, p) displayed in fig. [|, with normalization now adjusted according 
to the best-fit values given in the table. The diamonds indicate data from ZEUS 
JD and the squares data from HI ||. 

i) The curves show F2 computed from G as displayed in fig. (|i) using eqn. (|T6|); 
the dashed line is a linear fit to the HERA data; a) is the soft pomeron prediction 
(solid curves of fig. f|i)) and b) is the hard pomeron (dashed curves of fig. |]i)). 
The subasymptotic data from the NMC are also shown (crosses). The filled points 
correspond to data with p < 1.4. 

ii) The curves are now calculated from G as displayed in fig. |]ii), compared to 
the HERA data rescaled by Rp eq.fjTj]). Again a) is the soft pomeron prediction 
and b) is the hard pomeron. 
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